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ABSTRACT 

Previous observational studies of the infrared (IR)-radio relation out to high redshift employed any 
detectable star forming systems at a given redshift within the restricted area of cosmological survey 
fields. Consequently, the evolution inferred relies on a comparison between the average IR/radio 
properties of (?) very IR-luminous high-z sources and (ii) more heterogeneous low(er)-z samples 
that often lack the strongest IR emitters. In this report we consider populations of objects with 
comparable luminosities over the last 10 Gyr by taking advantage of deep IR (csp. Spitzer 24izm) and 
VLA 1.4 GHz observations of the COSMOS field. Consistent with recent model predictions, both Ultra 
Luminous Infrared Galaxies (ULIRGs) and galaxies on the bright end of the evolving IR luminosity 
function do not display any change in their average IR/radio ratios out to z ~ 2 when corrected for 
bias. Uncorrected data suggested ~0.3 dex of positive evolution. 

Subject headings: cosmology: observations — galaxies: active — galaxies: evolution — infrared: 
galaxies — radio continuum: galaxies — surveys 



1. INTRODUCTION 

The IR/radio properties of galaxies at successively 
higher redshift have been probed in the past decade 
using either st atistical samples from cosmological sur- 
vey fields (e.g. lAppleton et all 120041; IFraver et all 120061: 
ISargent et all 120101 hereafter: S10), the stacki ng tech- 
nique (e.g. ICarilli et al.ll2008HIvison et al.ll2010D or ded- 
icated sampl es of specific objects (e.g. sub-mm galax- 
ies (SMGs); IKovacs et all [2005 IHainline et all l200l 
iMichalowski et al.ll2009() .Evolutionary studies, all based 
on samples poorly matched in terms of bolometric lu- 
minosity at low and high redshift, have provided con- 
flicting results, conc luding that the local IR-radio rela- 
tion either does (e.g. IGarrett]|2002T: lAppleton et al.H2004l: 
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Ibar et al.ll2008t iGarn et al.ll2009{: S10) or does not (e.g. 
redshift. 

Recently, predictions have been made for the redshift 
evolution of the IR/radio properties of star-forming 
galaxies having different luminosities and geometries 
(e.g. compact starbursts and normal sta r-forming disks; 
iMurphvl l200l lLacki fc Thompson! 120091 ) . The current 
generation of IR and radio observatories can directly de- 
tect the brightest of these systems over a significant frac- 
tion of Hubble time, provided that sufficiently large cos- 
mological volumes are sampled. Here we make use of the 
VLA and Spitzer coverage of the 2 deg 2 COSMOS field to 
construct (cf. Ej2j) a volume-limited sample of ULIRGs at 
z < 2 that allows a direct comparison of observations and 
theory. Our findings are presented in SJ3] and discussed 
in §|2 

We adopt the WMAP-5 cosmology (fi m = 0.258, 

^A + ^m = l and i?o = 71.9 kms _1 Mpc _1 ; 

IDunklev etld1[2009t ). 

2. DATA AND SAMPLE SELECTION 
2.1. IR and Radio Measurements 
The 1.4 GHz map (jSchinnerer et al.l 120071) of the 
2 deg 2 COSMOS field reaches an average sensitivity of 
-0.017mJy/beam (FWHM = 2.5"). Here we use the 
VLA-COSMOS 'Joint' Catalog (Schinnerer et al. 2010, 
subm.) containing ^2900 sources detected with S/N > 
5. Spitzer/MIPS im aging by the S-COSMOS project 
(|Sanders et al.1 120071 ) achieves a resolution of 5.8" 
(18.6") and a la point source detection limit of ^0.018 
(1.7) mJ y at 24 (70) /xm (f or details see ILeFloc'h et all 
120091 and IFraver et alJl2009f ) . The depth of the 24 /zm ob- 
servations exceeds that at 70 /jm and 1.4 GHz by roughly 
a factor of seven in terms of equivalent IR luminosity (cf. 
Fig. 1 in S10). At an equal detection significance level 
(3ct), the 24 /im catalog consequently is roughly 20- fold 
larger than the 70 fim source list (~50,000 vs. 2,700). 
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Spitzer detections were matched to the VLA-COSMOS 
sources using a search radius of FWHM /3 for either IR 
filter. Ambiguous radio-IR associations were removed 
from the sample. Since the 1.4 GHz catalog is restricted 
to sources with S/N > 5, we re- analyzed the radio map 
at the position of unmatched IR sources and added all 
resulting detections with S/N > 3 (^2100 objects) to the 
sample. For more information on the band-merging and 
the flux distributions of the (matched and unmatched) 
IR and radio sources we refer to the detailed description 
in S10. 

We use the joint flux information at 24 and 70 /im to de- 
termine - given the known redshift (cf. § 12.21) - the best- 
fitting synthetic IR spectral energy distribution (SED) 
and thence the IR luminosity, Ltt r = £(8 — 1000 /xm) . 
As described in S10 (§4; see also Murphy et al.1 [20091 
for additio nal details on th e SEP fitting), templates ac- 
cording to iCharv fc Elba3 (|200l are used for galaxies 
directly detected in both MIPS filters. For sources only 
detected at 24 /jm we also fit the Spitzer p hotometry (in- 
cludin g the 70 /im upper flux limit) with iDale fc Heloul 
<|2002f ) template SEDs and define the best estimate of the 
IR luminosity as the average Ltir from the two separate 
fits. Inferring Ltir from only two bands at wavelengths 
shorter than the peak of the IR emission is expected to 
lead to uncertainties of a fact or two to five (jMurphy et al.l 
l2009HKartaltepe et al.ll2010h . Although our estimates of 
Lttr are thus less precise than those presented by, e.g., 
Ilvison et al.1 (|2010l ) for a similar study, there is no consen- 
sus in the literature that they are systematically biased 
to low or high fluxes (cf. discussion in S10, § 6.5). 
The IR/radio properties of our sources are charac- 
terized by the lo garithmic TIR/radio flux ratio Qtir 
(|Helou et al.lfl985l) : 

<?tir = bg ^ ?5 x 1Q12 w j - log ) . (1) 

The rest-frame 1.4 GHz luminosity Li.4GHz is 

ii^GHzfWHz- 1 ] = ^4r£ £,(1.4 GHz) , (2) 
(l + z) 1 a 

where 5^(1.4 GHz) is the observed integrated radio flux 
density of the source and Dl(z) the luminosity distance. 
The X-correction (1 + z) - ^ 1- ") depends on the spec- 
tral inde x of the synch rotron emission, which is set to 
a = 0.8 (jCondonlll992l ). Given the mean spectral slopes 
typically measure d for faint extra galactic radio sources 
(0-4 < a < 0.9; e.g. llbar et afll2"00ll our values of Li.4GHz 
should be accurate to within 40 (70)% at z~ 1 (2). The 
main contribution to uncertainties on <7tir thus stems 
from errors on £tir- 

2.2. Distances and Source Classification 

Optical data and photometric r edshif tj^l are taken 
from the catalog of lllbert et all (120091 ) . The wave- 
length range spanned by these observations (30 broad, 
medium and narrow band filters) extends from 1550 A to 

15 The photo-z dispersion a(&z/(l+z)) is 0.007, 0.013 and 
0.051 for sources at z < 1.25 with i A B < 22.5, iab G [22.5, 24] and 
«AB > 24, respectively. At higher redshifts the accuracy of the pho- 
tometric redshifts decreases by a factor of ~3. 



8 /im. iCapak et all ([20071 . 12008ft provide a complete de- 
scription of these observatio ns. Spectroscopic redshifts 
from the zCOSMOS survey (jLillv et al.l [2009]) or Mag- 
ellan /lMAf2S_a^ Keck/Dcimos follow-up observations 
(e.g. iTrump et al.l 120091 : Kartaltepe et al., in prep.) are 
available for ^25% of our sources, most of which lie at 
z < 1. (See values of f S pectro(-z), the spectroscopically ob- 
served sample fraction, in Fig. 2.) The quality of dis- 
tance measurements is assessed using spectroscopic con- 
fidence flags or the width of the photo- z probability dis- 
tribution in the case of photometric redshifts (see S10, 
Appendix). Sources that do not fulfill the reliability re- 
quirements are excluded from the subsequent analysis. 
The optical photometry and spectroscopy was matched 
to the radio and IR catalog entries using a search radius 
of 0.6" and 1" (reflecting the larger uncertainty on the 
centroids of IR sources), respectively. As done for the 
band-merging of the Spitzer and VLA data, ambiguous 
optical- IR/radio associations are discarded. 
In Fig. [T] we plot the IR luminosities of our sources 
function of their redshift. Based on the range 
of accessible IR luminosities at each redshift we de- 
fine two populations for later investigation: (i) ULIRGs 
with Ltir > 10 12 L©, and (ii) all objects populating the 
bri ght end of the TIR lu minosity function (LF) derived 
bv lMagnelli et al.l (|2009[ ). The bright end is defined as 

Lrm > 4iR e) (*), wh ere L%£*\z) cx (l+z) 3 6±0 - 4 repre- 
sents the break in a double power law parameterizatiorF^I 
for the TIR LFs. Both selection approaches lead to a 
volume-limited sample of either ULIRGs or 'IR-bright' 
galaxies spanning the range z < 2. 

We divide our sample into star forming (SF) galaxies 
and active galactic nuclei (AGN) using a modification 
of the rest-frame optica l color-based method developed 
by iSmolcic et all (|2008f ). Optical-to-near-IR SEP fits 
with the package ZEBRA (jFeldmann et al.ll2006l ) provide 
rest-frame (u — A")-colors that are translated into the 
probability of 'SF-hood', Pr(SF), for each object in our 
sample (cf. § 3 in S10). SF systems are all galaxies with 
Pr (SF) > 0.5 (or (u-K)ab < 2.36). Galaxies with redder 
colors are regarded as AGN hosts. 

3. RESULTS 

To constrain the evolution of average IR/radio ra- 
tios we compute the median, ((/tir), in different red- 
shift slices. The selection threshold for ULIRGs and the 
IR-bright population lies well above the faintest acces- 
sible IR luminosities in the redshift range < z < 2 (cf. 
Fig. [IJ. Our samples are a mixture of (i) sources di- 
rectly detected at IR and radio wavelengths, plus (ii) 
24 /im-detected sources with only a 3 ct upper radio flux 
limit from the 1.4 GHz rms image. The corresponding 
IR/radio ratios are either well-defined (within experi- 
mental uncertainties) or lower limits and, when com- 
bined, form a 'censored' sample that is best analyzed 
with the tools of survival analysis. In the present case of 
one-sided censoring, the cumulative distribution function 
(CPF) of measurements of <?tir can be derived with the 

16 Consistent with the meas urement of Magnclli ct al.l H2009T ). 
ILe Floc'h et al.l ( 120051') and lCaputi et al.l (120071) report an evolution 

of (l+z) 3 ' 2 - - 2 and (l+z) 3,5±0 - 4 , respectively, based on a double 
exponential fit to the TIR LF at z < 1. 



IR-bright galaxies on the IR-radio relation 
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iKaplan fc Meierl ([1958) product limit estimator. As it 
is normalized (i.e. runs from zero to unity), the median 
corresponds to that (/tir for which the CDF is equal to 
0.5. 

We derive (q T m) for our sample of 1,692 SF ULIRGs 
and for all 3,132 COSMOS ULIRGs (SF and AGN). 
Fig. [2a| shows the accordingly normalized CDFs. At 
<?tir ^ (<7tir), the CDFs approach a non-zero value re- 
flecting the number of lower limits on (?tir that may 
exceed the largest uncensored measurement in that bin. 
Due to the comparatively shallow VLA observations, IR 
sources without a directly detected 1.4 GHz counterpart 
become more frequent as redshift increases. Therefore, 
the width of the redshift bins chosen for the construc- 
tion of the CDFs compromises between a split of the 
studied redshift range z < 2 into regular intervals and 
the aim to sample th e di stribution function down to 
the median. In Fig. [2b] we show the CDFs for IR- 
bright COSMOS sources (3,004/5,657 in the SF/total 

sample, respectively) that satisfy Ltir _ L^^ c \z). The 
larger number of faint (cr) luminosity sources allowed us 
to divide the range z<1.6 into thinner slices than was 
done for ULIRGs. Moreover, there was a sufficiently 
large number of objects of this class even at low redshift 
(0<z<0.2), whereas the closest ULIRG in our sample 
lies at z-0.4 (cf. Fig. p. 

The values of (<7tir) reported in Table[T]increase at z 3> 1. 
In S10 (Fig. 17 and Table 6) we found that - although 
average IR/radio ratios display little evolution with red- 
shift - the dispersion <r qTIR in the COSMOS data in- 
creases fairly abruptly at z^l.4. This is primarily due 
to increased uncertainties in <7rrR, but might also hide 
a small intrinsic increase in the dispe rsion. The scatter 
fT 9TIR directly influences the shift (e.g. lKellermannlll964| ) 



A 9b 



ln(10)(/3-l> 



2 

<?TIR 



(3) 



between the average IR/radio ratio of flux-limited sam- 
ples selected at IR and radio wavelengths ((3 is the power 
law index of the differential source counts dN /ds oc S'-' 3 ). 
S10 showed that eq. ([3]) predicts the actual offsets 
present in the COSMOS data remarkably well. Because 
of the higher sensitivity of the 24 /j,m observations, the 
present sample is effectively IR-selected. Eq. © allows 
us to compensate for the relative offset between medians 
at high and low redshift that arises artificially due to the 
increased scatter in our data at z> 1.4. In doing so, we 
(i) use that a qT1Ii (z < 1.4) w 0.35 and ct 9tir (1.4 < z < 
2) « 0.75 (cf. S10, Table 6), and we (ii) assume that the 
observed flux densities of galaxies at z > 1.4 primarily lie 
in a range of sub-Euclidean source counts wh ere (3 w 1.5 
(e.g. ICharv et al] 120041: IPapovich etHI [200l . The re- 
sulting correction (to be subtracted from the medians at 
z > 1.4) is 

Dg corr .=Ag bias (1.4 < z < 2) - Aq bias {z < 1.4) (4) 
= ln(10) [(1.5 - 1) x 0.75 2 - (2.5 - 1) x 0.35 2 ] 
-0.22 , 

with an associated uncertainty (owing to the errors on 
(7 9TIR and (3) of approximately 0.13. Our step function- 
like correction neglects that sources at a given redshift 
may be drawn from a flux range with continuously vary- 
ing (3. This is the simplest possible form that allows us 



to correctly compensate for an apparent, spurious offset 
in ((2tir) between the limits of our investigated redshift 
range. 

In Fig. |3] we plot (<7tir) vs. redshift and relate these 
medians with the best-fitting evolutionary trend of the 
form (9tir(z))/( 9tir ( z= o)) oc (1 + z) 7 for each of our (sub- 
samples. To constrain the fit at low redshift we add a 
low-z data point (both for the ULIRG and the IR-bright 
sample) based o n the complete IRAS-selected sample of 
lYun et all (|2001f ). The values of (/fir given bv lYun et al.1 
(|2001l ) were converted to <7tir by boosting their IR flux 
by a factor of twcFl. the average diffe r ence between 
the mean (?tir and <7fir foun d bv [Bell (l200l bmce 
SF systems in the sample of lYun et all (|2001l ) could 
not be identified following the procedure employed for 
the COSMOS galaxies (see § 12. 2p . our reference sample 
for SF galaxies simply consists of all local sou r ces w ith 
£i.4GHz < 10 24 W/Hz, in keeping with lCondonl (fl989l ). 
The two panels of Fig. [3] show the evolutionary trends 
for SF systems (top) and all galaxies (SF and AGN; bot- 
tom). To account for the asymmetric error bars, we have 
drawn 1001 values from within the 95% confidence region 
of each mediarF^I and then fit the evolutionary trend 1001 
times based on random combinations (without replace- 
ment) of the resampled medians. The final best-fit values 
given below are the medians of the parameter distribu- 
tions thus obtained. 

If the uncorrected high-z medians arc included in the evo- 

+0.08 , 



lutionary fit, we find an exponent 7 = 0.091^/0.11 



I +0.06 
-0.05 

for ULIRGs (SF and total sample, respectively) and 
7 = 0.13 ± 0.06/0.12 ± 0.04 for the IR-bright galaxies. 
(Errors delimit the 95% confidence interval.) This would 
imply a doubling of the average TIR/radio flux ratio from 
redshift to 2, with most of the evolution happening 
at z>1.4. By taking the corrected medians at z>1.4 
into consideration, we find an evolution of the average 
IR/radio ratios of ULIRGs according to (1 + z )-0-Oi±o.06 
(identical for SF and all ULIRGs). Similarly, the trend 
for the IR-bright population is consistent with being zero, 
both for the total sample (7 = O.OSlg^) and the SF sub- 



set where 7 = 0.04 



+0.06 
-0.05- 



4. CONCLUSIONS 

We have presented the first investigation of the evolution 
of the IR-radio relation out to z ~ 2 for a statistically sig- 
nificant, volume-limited sample of IR-luminous galaxies. 
This advance became possible thanks to two factors: (i) 
the large area and deep mid-IR coverage of the COSMOS 
field, and (ii) the inclusion of flux limits in the analysis 
with appropriate statistical tools. 

At redshifts z < 2 the median TIR/radio ratio of ULIRGs 
remains unchanged if we compensate for biases. On 
the most basic level this implies that their magnetic 
fields, B, are sufficiently strong to ensure that cosmic ray 
electrons predominantly lose their energy through syn- 

17 Due to the well-con strained mean IR/radio r atios in the 
samples of lYun et all :|21MU ((^fir) —2.34+0.01) and lBeH JUS) 
((<?tir) = 2.64+0.02) this average correction factor is accurate to 
within a few percent. 

18 When upper confidence limits are 00, we set the upper error 
bar to twice the lower one. This eases the calculation of 7, while 
reflecting that upper and lower confidence interval are generally 
similar down to the 60 th percentile of the CDF. 
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chrotron radiation (rather than inverse Compton scat- 
tering off the CMB). Regarding the 0.3 dex increase of 
the uncorrected evolutionar y signal as an upper limit 
implies that B>30fiG (e.g. IMurphvl 120091 Fig. 5), as 
expected for compact and strong starbursts. This con- 
clusion applies to both SF systems and optically selected 
AGN hosts, consistent with the similar mean IR/radio 
ratios reported for these two classes of objects in S10. 
Our finding agrees with theoretical and numerical ex- 
pectations that ULIRGs shoul d follow the local IR-radio 
relation until at least z ~ 2 (jLacki &; Thompson! 120091 : 
IMurphvl 120091 ) . Moreover, it suggests that the lower 
IR/radio ratios frequently reported for high-z SMGs (e.g. 
Kovacs et all 120061; iValiante et al.l l2007t iMurphv et al. 
2 9 t I MichalowskTc t al. 2009], but see also lHainline et ah 
2009J) are not typical of distant ULIRGs in general. 
Our complete sample of 'IR-bright' galaxies - the pop- 
ulation that resides on the evolving bright end of the 
TIR LF - links high-z ULIRGs to normal IR-galaxies 
(log( L Tm/L ) > 10.5) in the local universe. The fact that 
the average IR/radio ratio of the latter is very similar to 
that of ULIRGs demonstrates that the similar IR/radio 
properties of existing SF samples at low and high rcdshift 
are not the fortuitous consequence of comparing objects 
in different luminosity ranges. While distant starbursts 
follow the same IR-radio relation as local sources, this has 
not yet been ascertained for galaxies with moderate SF 
rates (<lOM /yr) that cannot be directly detected with 
current ra dio and far-IR faciliti es. The recent stacking 
analyses of lSevmour et al.l (|2009f ) and llvison et al.1 (|20ld ) 
measured a steady decline of average IR/radio ratios that 
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TABLE 1 

Medians of the CDFs shown in Fig. 2. Bracketed values of (<ztir) are the medians in high redshift bins before correction 

according to eq. jsj. 



sample/galaxy type (z) (<3tir) N (z) (<2tir) N 



ULIRG sample 'IR-bright' sample 



local (all sources) 


n 099+ - 020 
u.uyy_ 004 


2 674+ ' 020 




51 


0.022±0.001 


9 091 +0.016 
z ' Dzi -0.013 




1,109 


COSMOS 










0.169+.0.014 


9 tr 9 o+0.054 
z.ozo_ 061 




34 


(all sources) 










0.346±0.008 
0.518+.0.008 


9 7n o+0.041 
z. <uo_ 040 

2 629+ ' 067 




231 
272 




0.659±0.053 


2 622+ 179 




26 


0.705±0.004 


9 5-1^+0.017 
z.o±o_ 047 




486 




0.943±0.022 
i 900-1-0 noQ 


2 439+° 124 




52 
213 


0.906±0.005 
1.108+.0.006 
1.280±0.007 


9 t-c9+0.053 
z.jjz_ Q 034 

9 R1 9+0.034 
z - Dlz -0 .086 

2 619+° 170 
Z.U13_ 056 




754 
637 
487 




1.466±0.004 


2.520l 042 


(2.740) 


715 


1.498+.0.005 


9 p. 7fi +0.239 


(2.796) 


1506 




1.628±0.003 


9 585+ - 245 
z.ooo_ 057 


(2.805) 


553 












1.785+.0.004 


9 71 

z - ' io -0.144 


(2.935) 


627 


1.784±0.004 


9 71 K + °° 

z - ' lo -0.146 


(2.935) 


546 




1.981±0.003 


2 651+ 0105 
z.oo±_ 104 


(2.871) 


1,185 


1.980±0.003 


2 651 + 0140 
Z.DOl_ 102 


(2.871) 


875 


local (SF sources) 


095+ - 015 


9 7HQ+0.055 
z. ( uo _o.050 




47 


0.022±0.001 


9 fi 99+0.017 
z.uzz_q 012 




1,101 


COSMOS 










0.169±0.015 


2.538+g-lS 

— U. llO 




19 


(SF sources) 










0.340±0.014 


9 71fi +0.075 
z. 1 1D_ .064 




101 




0.624±0.062 


9 799 +0.349 
z - ,z -0.279 




10 


0.528±0.010 
0.705±0.006 


2 649+ ' 086 
z.oty_ 117 

2 ^40+ ' 092 
Z.O4U_ 035 




142 
242 




0.940±0.030 


9 4QO+0.156 
i ™-0 .212 




27 


0.915±0.007 


2 560+ 053 

Z.ODU_ 037 




394 




1.210±0.012 


9 i 7+0. 132 
Z - D1 ' -0.082 




119 


1.107±0.008 
1.273+.0.009 


9 K 79+0.069 
-0.061 

2.709lS° 147 




327 
263 




1.467±0.004 


z.oot_ 03g 


(2.754) 


412 


1.501+.0.006 


9 fifi/1+0.258 
Z.OO1_ 066 


(2.576) 


932 




1.628±0.004 


9 c;7fi+0.254 
z -°' D -0 .080 


(2.796) 


348 












1.786±0.006 


9 71 5+°° 
z - ' io -0.131 


(2.935) 


321 


1.786±0.006 


9 71 r + 00 
z - ' ±o -0.132 


(2.935) 


268 




1.984±0.004 


2 672+°° 
Z - D,z -0.121 


(2.892) 


591 


1.982±0.005 


9 71 9+00 
z - ' ±z -0.160 


(2.932) 


403 



Note. — Columns are: sample description (1), median redshift (with 1 a error) of redshift slice in which CDF is constructed (2; 5), median 
TIR/radio ratio with upper and lower 95% confidence intervals (3; 6), number of objects used (4; 7). 
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Fig. l. — : Infrared luminosities Ltir of our 24 /im- and/or 1.4 GHz-detected sources (see legend in the upper left 
corner) as a function of redshift. The dotted/dashed line tracing the lower edge of the measurements mark the minimal 
detectable Ltir as predicted (based on 24 /j,m sensitivity) by template SEDs. Dark (light) grey symbols denote sources 
with spectroscopic (photometric) redshifts. The solid black line shows the evolution of the characteristic luminosity 

•^tir °f the TIR LF at z< 1.3 (jMagnelli et al.ll2009f ). Its extrapolation to z^2 is schematically indicated by the 
dash-dotted segment. Red and green lines delimit the luminosity range of the two samples defined in § 12.21 
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Fig. 2a. — : Cumulative distribution functions (CDFs) of ULIRG TIR/radio ratios, specifying the fraction of the 
population with a ^tir larger than a given value on the horizontal axis. The panels show every second redshift bin 
from Fig. [3] 

Hatched areas: 95% confidence interval. Grey: CDFs of SF objects. Black: all active galaxies. The intersection of 
the curves with the solid horizontal line defines the median (qtir) (vertical dashed lines). f spec tro specifies the sample 
fraction with spectroscopically measured redshifts; f ce nsored is the fraction for which 9tir is constrained by a lower 
limit. 
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Fig. 2b. — : As Fig. [2a] but for the IR-bright population. 
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Fig. 3. — : Redshift evolution of the median logarithmic TIR/radio ratio (qtir) for IR-bright galaxies (Ltir > ^tir > 
green symbols) and ULIRGs (red). In the upper panel we consider the subset of SF sources, extracted from the entire 
sample of active galaxies (bottom). Transparent symbols: estimates of (c/tir} prior to correction for selection biases 
(see §[3j). The best- fitting evolutionary trends to the corrected (uncorrected) measurements of (<?tir) are reported 
using stro ng (transp a rent) dashed lines. They have been additionally constrained (open stars) at low redshift by the 
sample of lYun et al.l (pOOl . 

Both ULIRGs and IR-bright galaxies have constant average IR/radio properties out to z ~ 2 when correcting for bias, 
otherwise ~0.3 dex of positive evolution is found. 



